The propagator for the step potential and delta function potential 
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We present a derivation of the propagator for a particle in the presence of the step and delta 

function potentials. These propagators are known, but we present a direct path integral derivation, 

i * ■ 

based on the path decomposition expansion and the Brownian motion definition of the path integral. 

The derivation exploits properties of the Catalan numbers, which enumerate certain classes of 

lattice paths. 

PACS numbers: 03.65.Yz, 03.65.Nk 
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I. INTRODUCTION 



Calculations involving step and delta function potentials occur in many branches of 
physics. Step potentials can be used to represent 'hard wall' boundary conditions, and 
are also involved in tunneling calculations. Imaginary step potentials are useful in modeling 
particle absorbtion, for instance in arrival time problems [l|, |2j. Delta function potentials 
can be used to model point interactions, especially in the low energy limit where the details 
of the process are largely independent of the form of the scattering potential [3j. In this 
paper we present a derivation of the propagators in the vicinity of these potentials by making 
use of the path decomposition expansion \4, la, ly] , and then by appealing to the Brownian 
motion definition of the path integral. 

We wish to calculate the following propagator 



0(xi,T|x o ,O) 



x(T)=xi 



T>x e 



iS 



x(0)=xo 



where 



S 



dt 



mx 



V(x 



(1) 



(2) 



'o v 2 
and the potential V(x) will either be a step potential V9(— x), or a delta function potential, 

aS(x). (Note that throughout this paper we set h — 1.) To calculate these propagators, 

the path decomposition expansion [4J, |5|, |6| will be employed as follows. A typical path from 

Xo > to x\ < may cross x = many times, but the set of paths may be partitioned 

by their first or last crossing times. We therefore split every path into three parts: (A) a 

restricted part that starts at xq and does not cross x — 0, but that ends on x = at time 

ti, (B) an unrestricted part from x = to x = that may cross x = many times and, (C) 

a further restricted part from x = to x\ that does not re-cross x = 0, Fig. (CD). In terms 

of the first crossing time ti, 



g(x!,T\xo,0) = — - / dti g(x 1 ,T\0,ti) -^-(x,ti\x ,Q) 



2m 



dx 



(3) 



x=0 



And in terms of the last crossing time £2, 

i r T Qg r 
0(xi,T|x o ,O) = -— / dt 2 -Q L (xi,T\x,t 2 



g(Q,t 2 \xo,Q). 



(4) 



x=0 



From the above it is possible to derive a third decomposition, in terms of the first and last 
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t=T , X 1 
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crossing times 
5f(xi,T|x ,0) 



4m 2 



FIG. 1: A typical path from xq to x\ 
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dx 



(5) 



When X\ and Xq are both positive, there is an additional contribution from paths which 
never cross x = and the path decomposition expansion is, 

dg r 



% f 
g(xi,T\x ,0) = g r (xi,T\x o ,0) - — / dh g(xi,T\Q, t x ) 

2m J l 



dx 



(x,h\x ,Q) 



(6) 



x=0 



In the above path decompositions g r are the restricted propagators for intervals (A) and 
(C), when the particle does not cross the origin. Since in the regions (A) and (C) the 
potential is a constant, V (or zero) the restricted propagators in these regions are equal 
to the restricted propagator for the free particle, weighted by a factor e~ v( - tl ^ t °\ Now the 
restricted free propagator, gf >r , is given by the method of images expression, 



gfr{xi,h\x ,t ) = 9(±x 1 )9(±x ) (g f (x 1 ,t 1 \x ,t ) - g f (-x 1 ,t 1 \x , t )) 



(7) 



where the 9 functions ensure paths start and end on the same side of x = 0. So the restricted 
propagators may therefore be easily computed by weighting the free propagators in Eq. (J7J) 
according to the region of interest. The problem of calculating the full propagator therefore 
reduces to that of calculating the partial propagator for the interval where the path crosses 
the origin g(0, t 2 \0, h)- We will show how these partial propagators may be derived using 



the Brownian motion definition of the path integral [7|, |8( . The full propagators may then 
be obtained with the help of Eqs. ([3]l5]) . 

nnnn 

For the step potential the full propagator has been derived in Refs. |9|, I Id. Ill [12J, but 



we shall derive the partial propagator, and then direct the reader to Ref. [10] for details 
of the use of the path decomposition expansion to recover the full propagator. The partial 
propagator we will derive is given by 



For the de 
from Ref. [13] 



, x ( m \V2 (l-e- iVT ) 

9(0 ,T|0,0) = ->(-) L_i. (8) 

ta function potential we will derive the full propagator and we quote the result 



POD 

g(x 1 ,T\x ,0) = gf(xi,T\xo,0)-a du e~ amu g f {\xi\ + |x | + u,T|0,0) (9) 

Jo 



where 



gf(x h T\x o ,0) = {^f) 1/2 e im(xi - X ° )2/2T ( 10 ) 



is the free propagator. 



II. THE BROWNIAN MOTION DEFINITION OF THE PROPAGATOR 



We begin with a review of some of the details o: 



' the Brownian motion approach to 
8J. The first step is to switch to 



7. 



computing propagators. For more details see Refs. 

working with the Euclidean propagator g by means of a Wick rotation and we specialise 

immediately to the case of xq = X\ = 0. That is, we wish to calculate 

px(T)=0 

g(0,T\0,0)= / Vxe~ SE (11) 

Jx(0)=0 

where Se is the Euclidean action given by 

S E = J Q T dt(^f + V(x)y (12) 

This propagator may be viewed as a conditional probability density for a random walk on 
the real line. The second step is then to make this integral over paths into a concrete object 
by defining it as the continuum limit of a discrete sum on a lattice. 

To establish conventions and demonstrate the basic ideas we compute the case of a free 
particle, following closely the treatment in Ref. |8(. We consider a rectangular lattice 



with spacing in the time direction of e, and spacing in the x direction of rj, and consider 
propagation for a time T = 2en, so we have 2n steps in our paths (the reason for this choice 
is that it simplifies a number of later expressions, and avoids clumsy factors of 1/2). The 
conditional probability u(0, T|0, 0) to start at (0,0) and end at (0,T) is given by the number 
of paths connecting the start and end points, divided by the total number of possible paths. 
The set of all possible paths is bounded by the extremal paths that take n steps to the 
left/right, followed by n steps to the right/left, see Fig. ([2]). Since a path must have the 
same number of steps to the right as to the left to end up back at x = we find, 



w(0,T|0,0) 



1 (2n 



(13) 



2 2n v n 

The Euclidean propagator ~g is then defined as the continuum limit of uj2r\ where we take 
e, rj — > 0, n — > oo, keeping e/rj 2 = m and T = 2en fixed. That is, 



^(0,T|0,0):= lim (2^-^(0, T|0,0) 



m \ V 2 
2^f) 



(14) 



which is the expected result for the Euclidean free propagator. 



Limit of allowed 
x paths 




FIG. 2: A typical path from x = to x — 0. The dashed line shows the area within 
which all paths must remain. 



III. THE STEP POTENTIAL 

The propagator along the edge of a step potential is given by Eqs. (ITT]) and (fT2l . with 
V(x) = V6(—x). We can write this as 

g{0,T\0,0)= Vxexpt- dt^^-jexpt-V dt9(-x)j (15) 

which is similar to the free particle case except that paths are weighted by a factor e~ Vr 

where 

T = dt6(-x) (16) 

Jo 

is the length of time spent in x < 0. In the lattice case the corresponding conditional 
probability w(0,T|0, 0) is given by a sum of paths, each weighted by a similar factor. This 
may be written as 

n 

u(0,T\0,0) = —J2n k e- 2k * v , (17) 

fc=0 

where n^ is the number of paths spending a time 2/ce < T in the region x < 0. Expresions 
for these n^ are known and are in fact independent of k [141 ] . They are equal to the Catalan 
number 

n + 1 \ n J 
where 2n is the total number of time steps. We can see this as follows. First, note that 
the number of paths that never enter x < is given by C n , this being one definition of 
the Catalan numbers. Next consider the following mapping on any path spending a time 
2ek < T in x < 0, Fig.©. 

1. Start from t=0, and follow the path until it first crosses x — (if it doesn't cross then 
stop, the path is in the set of non-crossing paths.) 

2. Follow the path until it comes back to x = again, note the step at which this happens. 

3. Swap the section of path before this step with the section after it. 

4. The new path will now spend 2 fewer time steps in x < 0. 

By repeated application of this mapping, any path can be transformed into one which never 
crosses x = 0. The important point about this mapping however, is that it is bijective |l4j, 
which proves that the number of paths spending time Ike in x < must equal the number 




FIG. 3: Application of the rule to the path on the left produces the path on the right, 
with 2 timesteps less spent in x < 

of paths that never cross, for any value of k. This shows that n^ = C n for all k < n. So we 
now have 



^(0,^10,0) = ^J2a 



„ 2keV ^n sr^ 



-2keV 



k=0 

2^T 



2 2r 



fc=0 



(19) 
(20) 



1 - e-w 

This coincides with the free particle case if V = 0. Since we plan to take the continuum 
limit we can Taylor expand the exponential in the denominator to first order in e, and use 
the following useful assymptotic form for C n |l4( 



a 



4" 



to get, 



u(0,T|0,0) 



1 _ e -2e(n+l)V 



v / ?2en 3 / 2 V V 

Now take the continuum limit as in Section 2 to obtain, after some simple algebra, 

m\i/2(l-e- FT ) 



9(0,T|0,0)=(- 



VT 3 / 2 



which implies 



<?(0,T|0,0) 



/ m \i/2 (1 -e~ lVT ) 
Vhri) VT 3 / 2 

7 



(21) 



(22) 



This is our first result, the propagator along the edge of a step potential. The full propagator 



from xq to x\ may now be obtained by making use of Eq. (jSj) 10|. 



IV. THE DELTA FUNCTION POTENTIAL 

We now wish to evaluate Eqs. fill I) and (TT2T) with the potential given by V(x) = aS(x). 
In a similar way as for the step potential, we can write this as, 

#(0,T|0,0) = / Vxexp I- / dt^-Jexp(-a/ dt8(x(t)U (23) 

which is again similar to the free particle case except that paths are weighted by a factor 

e~ au where 

a= f dtS{x{t)) (24) 

Jo 

is the number of times a given path crosses x = 0. We model the delta function as a square 
potential of width rj, and height 1/rj so that every crossing of x = is weighted by a factor of 



e -ae/»?_ gj nce e j r ji — m ^ we can rewrite this as e~ r]am 15]. (a is usually positive, so that this 
represents suppression, but there is no need to assume this for the derivation.) We take the 
potential to be located on the left of x = so that to start and end at x = involves an even 
number of crossings. Such choices are only made for convenience and have no significance in 
the continuum limit. The conditional probability density u may be partitioned in a similar 
way to that for the step potential. However now the partitioning is with respect to the 
number of times a path crosses the square potential. Since the number of crossings will 
always be even we partition into classes of paths that cross 2/ < 2n times, so the conditional 
probability is 



1 
M (0,T|0,0) = — ^J(n,/)e" 2 ^ (25) 

1=0 

Where J(n, I) is the number of paths of 2n steps that cross the delta potential 21 times. It is 
known that these J(n, I) are given by the 21 th convolution of the Catalan numbers [16(, this 
can be demonstrated by writing a general path in terms of sums over non-crossing paths. 
These convolutions form the diagonal elements in Catalans triangle ly, [l7|, [l8| . We need to 
know the n th element in the 21 th diagonal from the right, which will give us J(n, I). 



From the formula for the elements of Catalan's triangle 161 ]. 



C ^ k)= fc!((n + l)! (26) 



it follows that 

J(nl)=c(n + ln-l)= ggM + JJ _ (2n\ n!n!(2Z + l) 

where we have extracted the binomial factor for later convenience. We cannot perform the 
summation in Eg. (1251) directly as we did for the step function potential, so we take the 
continuum limit first to leave ourselves with an integral. In order to do this we need a 
simpler form for the J's. We need the asymptotic form for n —>■ oo, but plotting Eq. (!27l) 
as a function of / shows that the dominant contribution comes from taking I —> oo as well. 
It is possible to derive the following asymptotic form, 

2n\2l _ f/n 



J(n,l)~( ~e-"», (28) 

\n J n 

and as this form shows, the maximum value of J occurs when I is of order ^/n. This is 

a consequence of the Brownian motion origin of the paths, which therefore have Hausdorff 

dimension 2 [19]. The statement that a typical path crosses the origin an infinite number of 

times is a consequence of this fractal nature of a typical path. 

If we use the asymptotic form for the J's, Eq. (|28|) and make the change of variable in 

the partitioning, Eq. ( |25l) . u = 2i]l we can turn the summation into an integral, 

1 /2rA r 2rjn u u 2 

m(0,T|0,0) = — / du -— e -^e- amu . (29) 

v ' ' ' ; 2 2n V n J J 2r] 2 n y ' 

We now take the continuum limit as in Section 2 to obtain 

/ m 3 \ 1/2 /-oo 2 

S(0,r|0,0) = i^-^j J duue-^e- amu 

/"OO 

= g f (Q,T\Q,0)-a du e- amu g f (u,T\0,0) (30) 

Jo 

which, after rotating back to real time, yields 

/"OO 

g(0,T\0,0) = g f (0,T\0,0) - a / du e- amu g f (u, T\0, 0). (31) 

Jo 

We can now use the path decompsition expansion to obtain the propagator for paths which 
start at xq and end at x\. There are 4 cases, depending on the signs of xq and x\, we shall 
present one case here, the others follow in a very similar fashion. First note that the free 
propagator obeys the following, 

g f (xi,T\x o ,0) = — I dti g f (x 1 ,T\0,t 1 ) — ^(x,ti|x ,0) 



(32) 



x=0 



il dh ^w^ T \ x ^ 



(7/(0,^1x0,0) (33) 

x=0 



provided sign(xo) = sign(-Xx). (gf r denotes the restricted free propagator.) (If the start 
and end points have the same sign then there are paths between the two that never cross 
x = 0, so the expression for the path decomposition expansion has to be modified as in 



[5 



Eq.® [5|.) 

Secondly note the following identities, 

g f (a,t 2 \b,ti) = g f (-a,t 2 \ -Mi) = 9f(0,t 2 \b - a,t x ) = g f (0,t 2 \a-b,ti). (34) 

which are just expressions of the symmetry of the free propagator. 

Consider the case where xq < and x\ > 0, noting that u > we first use Eq. ( 1321) to 
attach the leg from (x , 0), 

/■oo 

g(0,t 2 \x ,0) = g f {0,t 2 \x ,0) - a / due- amu g f {u, t 2 \x , 0) 

Jo 

= g f (0,t 2 \x ,0) - a / due" amu gf(0,t 2 \x -u,0) (35) 

Jo 

where we have used Eq. (|34"|) to obtain the second line. Since Xq < we have that Xq — u < 0, 

so we can attach a leg to (xi, T) using Eq. (1331) . 

5f(xi,T|x ,0) = g f (xi,T\x ,0) -a / due' amu g f (xi,T\x - u, 0) 

= g f {x 1 ,T\x ,0)-a / due-^^/dxil + |x | + m,T|0,0) (36) 

as expected. Calculation in the other cases proceeds in a very similar manner, and confirms 
that Eq. (13"6"j) is valid in all cases. 
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